Nuclear moments for the neutrinoless double beta decay II 



C Barbero^, F. Krmpotic^, A. Mariano^ 

Departamento de Fisica, Facultad de Ciencias 
Q\ • Universidad Nacional de La Plata, C. C. 67, 1900 La Plata, Argentina. 



^ . 

(N 



(N 
> 

O 

o 

O 






and 
D. Tadic 

Physics Department, University of Zagreb 
Bijenicka c. 32 -10.000- Zagreb, Croatia. 



Abstract 



The recently developed formalism for the evaluation of nuclear form factors in neutrinoless 
^ ; double beta decay is applied to ^^Ca, ^^Ge, ^^Se, ^^^Mo, ^^^Te and ^^Te nuclei. Explicit ana- 
lytical expressions that follows from this theoretical development, in the single mode model for 
the decay of ^^Ca, have been worked out. They are useful both for testing the full numerical 
calculations, and for analytically checking the consistency with other formalisms. Large config- 
uration space calculations are compared with previous studies, where alternative formulations 
^ I were used. Yet, besides using the G-matrix as residual interaction, we here use a simple 5-force. 
j^ ■ Attention is paid to the connected effects of the short range nuclear correlations and the finite 
nucleon size. Constraints on lepton number violating terms in the weak Hamiltonian (effective 
neutrino Majorana mass and effective right-handed current coupling strengths) are deduced. 
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1 Introduction 

During the last years we have developed a new formulation for the neutrinoless double beta 
iPPou) decay, based on the Fourier-Bessel multipole expansion of the hadronic current, and on 
the angular momentum recoupling. First, we did it for the mass term within the single mode 
model (SMM) |]l|. Later on, full QRPA calculations were done for this term [^. More recently, 
the same procedure has been applied to the evaluation of the so called "recoil term" in the 
charged Majoron emission [Q]. Finally, the complete formalism, including the right-handed {V + 
A) hadronic current, was presented 0]. The physical substratum in this development is the same 
as in the previous works on the same issue @, §, ^ ^, 0, namely, the same weak Hamiltonian 
was used. Thus, one cannot expect to get sensibly different results for the corresponding 
observables. Yet, we have succeeded in expressing all nuclear PPou moments in terms of the 
matrix elements of only three well-known one-body spherical tensor operators: 



'-n)i 



Yajm(^) = E^n+OA(A;r„)ljM(r. 

n 

SALJM(fc) = E^>nJA(A;r„)[o-„®FL(r„)]jM, 

n 
PLJM{k) = Yj Tn3L{krn) [Pn ® >L(rn)] JM, 



which have been around in nuclear physics for more than 40 years pO| , |TT|. This makes our 
formulation to be specially suitable for the nuclear structure calculations, and more simple than 
other formulations §, |, 0, |, i- 

In fact, the Fourier-Bessel multipole expansion has also been used by Vergados [§ and by 
Suhonen, Khadkikar and Faessler [^, as the staring point. However, the final outcomes for the 
nuclear matrix elements in these two theoretical developments are quite dissimilar, not only to 
our formulas, but also to each other. They are also different to the formulas derived by Haxton 
and Stephenson |^ and by Doi, Kotani and Takasugi @] within the closure approximation. As 
a consequence, the alternative formalisms cannot be confronted analytically, and the only way 
to test the consistency among them is by way of numerical procedures. 

As far as we know, numerical calculations, using the formalism from ref. 0], have so far 
been performed only for the neutrino mass term p|, |l^. Thus, to complete our study on the 
(3(3oi, matrix elements, in this paper we carry out calculations for several (3(3 decaying nuclei 
that are attractive from the experimental point of view {^^Ca, "^^Ge, ^'^Se, ^^^Mo ^^^Te and 



^^^Te). A comparison with similar studies is also made. 

As most of the previous studies were performed in the framework of the QRPA model 



121 O, O, IT^, |T6[, we will use here mostly the same nuclear structure approach. Only the 



^^Ca ^^^ Ti-decay will be discussed in a simple shell model, in order to compare our formalism 
with that of Vergados |Q]. Also, for the sake of comparison, throughout this work the bare axial 
vector coupling constant Qa = 1.254 will be used, although the effective value g"/^ = 1 is 



preferable in nuclear physics [17, 18 



The outline of this paper is as follows: 

In Section 2 we summarize the main results for the new formalism developed recently 0]. 

In Section 3 the decay of "^^Ca is discussed within the SMM and a comparison is made with 
the shell model calculation performed by Pantis and Vergados [|19|. Here we also derive the 
analytical expressions for nuclear moments, which we later use to test the numerical results. 
The competition between the effects of the finite nucleon size (FNS) and the two-nucleon short 
range correlations (SRC) is discussed as well. 

Section 4 deals with full QRPA calculations for ^^Ca, '^^Ge, ^^Se, ^^^Mo, ^'^^Te and ^^o^e. 
We compare them with two previous QRPA evaluations, namely with: 1) the results obtained 
by Muto, Bender and Klapdor (MBK) |jl^ , where the formalism of Doi, Kotani and Takasugi 
has been utilized, and 2) the calculation performed by Pantis, Simkovic, Vergados and Faessler 
(PSVF) |T^ in the framework of the formalism developed by Vergados ||^. In this section 
we also show the limits on the (3(3ou coupling constants that we deduce from the most recent 
experimental data [0, |^, |2|, ^, ^. 

Concluding remarks are pointed out in Section 5. 



2 Neutrinoless double beta decay formalism 

The PPoi, half-life is expressed in the standard form |]^: 

[To.(0+^0+ri = {m,)'C,n^+{\)'Cxx + {v)'C,, (2) 

+ {m^){\)CmX + {fn^){r])Cm'n + W{'n)Cxri, 

where (171,^) is the effective neutrino mass and (A) and (rj) are the effective coupling constants 
of the {V + A) hadronic currents. The coefficients 

Cmm = {Mp — Mgt) Gl, 

Cxx = Ml_g2 + ^Ml^g^-'^M2-M^+g^, 

Cr,r, = Ml^G2 + ]:Ml_g^ - lM2+M,_gs + M^gg + M^M^g-r + Mjgs, 
y y 

Cmx = {Mp-MaT)[M2-g3-Mi+g4, (3) 

Cmr, = -{Mp-MaT)[M2+g3-Mi_g4 + M^ge + Mpg5], 

Cx, = -2M2_M2+^2 + I [M2_Mi_ + M2+M1+] ^3 - '^M,_M,+g^, 

contain the combinations of the matrix elements 

Mi± = Mgt' -GMj.±3Mp>, 

M2± = MaTu.±Mp^--Mi^. (4) 

The kinematical factors ^^ are given by eq. (3.5.17) in ref. and the nuclear matrix elements 
read [§: 

Mp = (^) J2 {-ymoj{pn)Wjoj{p'n')n'fj{pnp'n';u;j.Jff'^{pnp'n'-,J:), (5) 

\9a/ J^pnp'n' 

Mar = - E {-lfy\^Lij{pn)WLij{pn)nfj^{pnp'n'-ujj.)f^''{pnp'n'-Jl), (6) 

LJI^pup'ri' 

Mp, = -2(^-] E ^''^'^\Jl\L){Jl\L)WJOJ{pn)WJOJ{p'n') 

\9a/ LJ^pnp'n' 

X nyLipnp'n'^uj.Jff^pnp'n'-^J:), (7) 



Mar' = 2 Y. i''^''^\L'l\L){L'l\L)yVL'ij{pn)yVL'ij{p'n') 

LL'J^pnp'n' 

X nY,L{pnp'n; uj j^) ff^ {pnp n ] J^), 



^- = 9Ar ^ z'^'-'^WLij{pn)ff''{pnp'n';J: 



X -/, 



6ll' - {Jl\L){Jl\L') 



WL'i.j{pn)nLL'{pnp'n'; uj. 



- 2V6LWiLJll; IL'){L1\L') J2 ^^i%ipn)nll,ipnp'n'; uj.) , (9) 



K=± 



Mr = 10 Yl i^^^'~^^L'^{lL\J'){lL\L')W{12LJ';lL')W{12JJ';lL') 

LL'J'J^pnp'n' 

X WL'iApn)Wj>u{p'n')nl]L{pnp'n';ujj.)ff\pnp'n'; J^), (10) 

"^ LJ^pnp'n' 

X Wjoj{pn)W.nj{pn)n]j\{pnp'n';uJj.)ff\pnp'n'; J^). (11) 

Here R is the nuclear radius, M^ is the nucleon mass, the index a labels different intermediate 



states with the same spin J and parity ir, J = v2J+T and {L1\J) is a short notation for the 
Clebsh-Gordon coefficient (L010| JO). The angular momentum coefficients are: Q 

L S J y , 

l-n 2 3n ) 
X W{y^ln3n\J)W{LJlnTlln]llp). (12) 

and the two-body radial integrals are defined as 

■RTUpnp'n'-Uj.) = R J dke+^v{k;u;j.)Rl{pn-k)Ri,{p'n'-k), k' = 0,1, ±, (13) 

where 

v{k-..j,) = l^^^^y -J, = Ej, -\{E^ + E,) , (14) 



"'^We use here the angular monicntuni coupUng \{^,l)j)- 



is the "neutrino potential", and 

R'l{pn;k) = Rl{lp,npjn,nn;k) = Unp,i^{r)un^,i„{r)JL{kr)r'^^'''dr, 

d , 21^ + 1 ±1 



?(±) 



R'^>{pn;k) 



""-'^^'Hrf^^ 



2r 



Un,,,iAr)3Liqr)r dr. 



(15) 



There are also two additional matrix elements, namely Mp^ and Mgtui, which are obtained 
from Mp and Mq^ by the replacement 



v{k;ujjn) ^v^{k;ujjn) = — 



1 



7i{k + uj.y 

Finally the two-body state dependent particle-hole (ph) density matrix 



(16) 



(17) 



contains information on the wave functions of the initial (|0^)), final (|0t)), and virtual inter- 
mediate {\Ja)) states. 

In particular, within the QRPA formulation, and after solving both the BCS and the RPA 



equations for the intermediate (iV — 1, Z + 1) nucleus ||25|, the two-body density matrix becomes 
ff^{pnp'n'; J^) = \unVpXj.{pn) + MpW„lj^(pn)] \up'Vn'Xj.{p'n') + Un'Vp>Yj^{p'n')\ , (18) 
where the notation has the standard meaning |^, |2S . 



3 ^^Ca ^^^ Ti decay within the single mode model 

In the SMM the virtual states in the intermediate nucleus ^^Sc are: [0/7/2(^)0/7/2 (^)]j+, where 
J"*" = 0"^ ■ ■ ■ 7^ 1^ . When the harmonic oscillator radial wave functions are used and the 
excitation energy uj^r is taken to be zero, we can go a step further in the analytical calculations. 



Table 1: Radial integrals TZI^j {pnpn; cjjn) for the /3/3oi/ decay in ^^Ca. The excitation energy 
ujjn is taken to be zero, and harmonic oscillator radial wave functions were employed with the 
oscillator parameter u = Mcu/h. Short notation C = ^^p°°^^^ has been used. 

L J CvTZZ CuTZy^ CTZfj CTZ'jl CTZ^jl 



37230 12870 

18615 9805 -18615 
-690 
-690 
11835 16585 -11835 
4734 12870 

11835 7065 -11835 
6030 
6030 
8415 10485 -8415 
1870 12870 

8415 5445 -8415 
8910 
8910 
6435 7425 -6435 
990 12870 

6435 -6435 



The results for the radial integrals TZ1^,{pnpn; ujj£) are listed in Table |l|, and after perform- 
ing the summations on L, L' and J', as indicated in eqs. (|])-(0), we get: 
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Mgt' 
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\j' 



(19) 



where ff^{J^) = ff^{pnpn; J+). The coefficients Ax{J^) are given in Table 0. It should be 
noted that: i) in the SMM the matrix element Mp is always null, independently of the value 
for the excitation energy tujj, and ii) while the Fermi matrix elements arise only from even 
multipoles, the remaining matrix elements only come from odd multipoles. 



Table 2: Coefficients 


Ax{J^) for the matrix elements 


given 


by 


eq. 


m- 


J+ 980A^(J+) 


63700Agt(J+) 


245Ai 


niJ"-) 


245A2«(J+) 


47775At(J+) 


0+ 8687 
























1+ 




-746499 




-738 


- 


-1218 






54327 


2+ 1315 
























3+ 




-117247 




-552 




-720 






36751 


4+ 459 
























5+ 




-72575 




-870 




-510 






21275 


6+ 175 
























7+ 




-91875 




-2450 











18375 



In the shell model (SM) the virtual states | J+) in the Sc nucleus and the wave function 
|0t) in the final nucleus ^^Ti, are described, respectively, as the one-particle one- hole and 
two-particle two- hole excitations on the ground state |0^) in '^^Ca. That is: 



l^+) = \[f7/2ip)fr-Mn)]J^), J+ = 0M+.-.7+, 

|o;) = Ea^lO/)l[/r/2(p)]^^[/772H]^+;o+), i' 



0+,2+,4+,6~ 



(20) 



and the two-body density reads [^] 



Table 3: Amplitudes (Ot| J"*") of the ground state wave function in ^^Ti and the corresponding 
coefficients p^^i^J^). 



J+ (0||J+) pP''{J+) 



0+ 


0.9433 


-0.0005 


1+ 


— 


0.0426 


2+ 


-0.3126 


0.1159 


3+ 


— 


0.2098 


4+ 


-0.1092 


0.3072 


5+ 


— 


0.3533 


6+ 


0.0231 


0.3017 


7+ 


— 


0.1186 



To get close to the results obtained by Pantis and Vergados [^ as much as possible, we use 
the same wave function for the state |0t) that was employed in their work. This wave function is 
listed in Table § together with the resulting values for the densities if^{J^). Note that p^'^(O^) 
and //^(l^) are relatively small because of the restoration of isospin and SU(4) symmetries, 
respectively. Their single-particle (Hartree-Fock) values, obtained from (Ot|0^) = 1, are far 
larger than the values shown in Table ^ (p^'^(0+) = — //^(1+) = —0.25). On the other hand, 
they are identically null when these symmetries are totally restored. 

The short-range correlations (SRC) between the two nucleons are taken into account via 
^ It is worth remembering g that, in the SMM, eq. ( [18|) reduces to 



o^\j^)=p'^csiJ^ 



u:o + G{J+) 



where f^Bcs^'-^^) ^ UpVnUnVp is the BCS value for the two-body density, ojq is the unperturbed proton- neutron 
quasiparticle energy, ojj are the QRPA energies, and G{J~^) = G{pnpn; J^) is the particle-particle matrix 
element. Thus we see that, within the SMM, the last factor in this equation plays the role of the effective charge 
for the PPov decay, induced by the QRPA correlations. 



the correlation function |2^ 

fsRc{r) = l-Ukcr), (22) 

where kc = 3.93 fm"^ is roughly the Compton wavelength of the w-meson. The finite nucleon 
size (FNS) effects are introduced in the usual way, i.e., by the dipole form factors in momentum 
space: 

,2 \2 



{gv,A)FNs — gv,A 



A^ 



A2 + P 



(23) 



with A = 850 MeV. The corresponding modifications of the neutrino potentials are shown in 
refs. i,i. 



Table 4: Nuclear matrix elements for the decay ^^Ca -^^^ Ti within the single-mode shell- 
model calculations. We have used uj^ = and four different results are presented: 1) (bare) 
no correlations and no nucleon form factor, 2) (FNS) no correlations but with nucleon form 
factor, 3) (SRC) short range correlations but without nucleon form factor, and 4) {FNS+SRC) 
correlations and nucleon form factor. 



M, 



GT 



Mp M 



GTUJ 



Mp 



M 



gt' 



Mp/ 



M„ 



Mr 



Present Results 

bare 

SRC 

FNS 

FNS-FSRC 



-1.168 0.177 -1.168 0.177 -1.168 0.177 -1.435 0.330 

-1.080 0.159 -1.080 0.159 -0.657 0.073 -0.105 0.284 

-0.960 0.134 -0.960 0.134 -0.644 0.066 -0.929 0.312 

-0.947 0.130 -0.947 0.130 -0.574 0.052 -0.796 0.309 
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bare 
SRC 

FNS 
FNS+SRC 



-1.216 0.185 -1.216 0.185 -1.216 0.185 -2.178 0.344 

-0.859 0.108 -0.856 0.108 -0.841 0.105 -0.115 0.346 

-0.986 0.134 -0.986 0.136 -0.635 0.063 -1.344 0.322 

-0.731 0.117 -0.731 0.098 -0.532 0.055 -0.324 0.330 



Present results are confronted with those obtained by Pantis and Vergados [|T^] in Table |^. 
They used a somewhat different approximation for the SRC and therefore it is plausible that 
our matrix elements do not fully agree with theirs in the second and fourth case. In the other 
two cases, they should be identical, but they are not! The difference is particularly pronounced 
for the recoil matrix elements M^. The reason for the discrepancies could be the values used 



for the harmonic oscillator parameter u = Mu/h and the nuclear radius R; we have utilized 
z/ = 0.916A-1/3 fni-2 and R = 1.2A^/^ fm. Q 

Anyhow it is worth noting that in both calculations the FNS effects and the SRC act 
coherently on the Fermi (F) and Gamow- Teller (GT) moments, in the sense that their combined 
effects always diminish them more than when they are acting individually. This, however, does 
not happen with M^, in which case the FNS+SRC values turn out to be significantly larger 
than the SRC ones. The explanation for this somewhat curious behavior of the recoil matrix 
element was given by Tomoda et al., [^ ^ and is as follows. The contribution of the weak 
magnetism in (H) can be decomposed into the central and tensor parts 0. The central part is 
the dominant one, and within the closure approximation and for cujir = 0, it can be rewritten 
in the form: ^ 

MUhare) = -^^^^(^| E^^^n^'^™ " ^'nK^^ - r.)k), (24) 



This matrix element is totally killed by the SRC p^ ) and therefore 



M^aiSRC) = -^J^{F\J2r^r^a^ ■ aJiv^ - r„)/5«c.(r„ - r„)|/) = 0. (25) 

OMjv 5' -4 mn 



The k'^ dependence of the form factors (p3D distributes the 5-function over a finite region |j^, 26| 
i.e., 



e/ \ FATS Ji _Ar 

o47r 



1 - Ar + -{Arf 



(26) 



Consequently, the matrix element (p^) decreases {Mj^dFNS) < Mjidbare)) and Mj^ciFNS + 
SRC) ^ 0. 



•^Our definitions for the nuclear matrix elements Mp,Mp/,Mp^ and Mj, agree with those of Pantis and 
Vergados Q only for g^ = 5v As we used here Qa = 1.254, the results hsted in their Table 1 have been 
renormalized accordingly. 

''The following relation has been used: 

S{r-r')= = ^J^^UmLir') f k'dkjiikr)jiikr'). 



10 



4 QRPA calculations 

We have employed a residual 6-ioYce V - 

Vs and Vt for the particle-hole, particle-particle and pairing channels 



A7f{vsPs+VtPt)S{r), with different strength constants 

2|, H. The single- 



pair I 



particle energies, as well as the pairing parameters v^°'^^{p) and v 
fitting the experimental pairing gaps to a Wood-Saxon potential well. 



n) 



have been fixed by 



Table 5: QRPA results for the nuclear matrix elements that include both the FNS and SRC 
effects. An average excitation energy (ujjn) of 5.0 MeV has been used in the present evaluation. 



Nucleus 




Mgt Mp 


Mgtu 


Mp^ 


Mgt' Mpi 


M« 


Mr 


Mp 


^^Ca 




















present 


-0.953 0.376 


-1.010 


0.361 


-0.022 0.203 


-1.888 


-0.033 


0.075 




ref. |15| 


-0.785 0.367 


-0.830 


0.343 


-0.765 0.395 


-1.522 


0.166 


-0.131 



76 



Ge 



82 



^e 



100 



Mo 



128 



Te 



130 



Te 



present 


-2.845 


0.749 


-2.864 


0.723 


-0.837 


0.371 


-4.863 


-0.065 


-0.889 


ref. 14 


-3.014 


1.173 


-2.912 


1.025 


-1.945 


1.058 


-3.594 


0.612 


0.530 


ref. ig 


-2.929 


0.111 


-2.683 


0.111 


-3.154 


0.102 


-7.423 


0.714 


-3.360 



present 


-2.717 


0.800 


-2.769 


0.771 


-0.603 


0.398 


-5.147 


-0.061 


-0.754 


ref. 14 


-2.847 


1.071 


-2.744 


0.939 


-1.886 


0.966 


-3.343 


0.789 


0.500 


ref. |15| 


-2.212 


0.018 


-2.124 


0.029 


-2.323 


0.009 


-3.700 


-0.175 


0.108 



present 


-2.155 


0.972 


-2.363 


0.935 


0.354 0.493 


-6.150 


-0.233 


1.265 


ref. 14 


-0.763 


1.356 


-1.330 


1.218 


1.145 1.161 


-4.528 


0.823 


-1.182 


ref. |15| 


-0.615 


0.471 


-0.420 


0.436 


-0.722 0.512 


-0.930 


0.293 


2.393 



present 


-3.417 


1.019 


-3.476 


0.980 


-0.835 


0.451 


-6.354 


-0.136 


-0.560 


ref. ig 


-3.103 


1.184 


-3.011 


1.047 


-1.999 


1.054 


-4.371 


0.583 


0.483 


ref. |15| 


-2.437 


0.044 


-2.179 


0.029 


-2.673 


0.054 


-1.522 


0.748 


-3.412 



present 


-3.225 


0.978 


-3.271 


0.938 


-0.819 


0.448 


-5.934 


-0.118 


-0.560 


ref. ig 


-2.493 


0.977 


-2.442 


0.867 


-1.526 


0.860 


-3.736 


0.574 


0.387 


ref. |15| 


-2.327 


0.009 


-2.083 


-0.002 


-2.553 


0.016 


-5.445 


0.656 


-3.376 
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As already mentioned, the proton and neutron gap equations have been solved for the 
intermediate (A^ — 1, Z + 1) nucleus as in ref. |^, and we deal only with one QRPA equation. 
Note that in this procedure we avoid the problem of overlapping of two sets of the same 
intermediate states generated from initial and final nuclei |]T^ . 



Table 6: The coefficients that appearing in eq. (0) (in units of yr ^), evaluated with the matrix 
elements given in Table ^. We have used the kinematical factors from ref. [§ . 



Nucleus 




^rnm ^XX 


C 


CmX 


^ rnri 


Ca»7 


^^Ca 
















present 


1.13 10-^3 7.42 10- 


-^^ 1.47 10"* 


-9.94 10-1^ 


2.56 10-11 


-8.10 10-13 




ref. m 


1.07 10-13 3 gg 10- 


-13 6.63 10-10 


-4.75 10-1^ 


-5.20 10-12 


-3.43 10-13 


7^n. 















present 8.27 lO'i^ 1.26 10-i3 8.20 10-^ 
ref. 1.12 10-13 1.36 10-i3 4.44 10-^ 
ref. [Ill 7.33 lO-i^ 1.12 10-i3 3.22 10"^ 



-4.61 10-1^ 
-4.11 10-1^ 
-4.49 10-1^ 



2.56 10-11 

2.19 10-11 

-1.54 10-11 



-1.57 10-13 
-4.99 10-1^ 
-2.11 10-13 



82 



^e 



present 



ref. [0 
ref. [Illl 



3.48 10-13 1.14 10-12 3.65 10-^ 
4.33 10-13 1.01 10-12 1.54 10-s 
1.75 10-13 4.78 10-13 1.53 10-^ 



-2.47 10-13 8.83 10-11 
-1.60 10-13 6.37 10-11 
-8.77 10-1^ -1.31 10-11 



-1.39 10-12 
-3.84 10-13 
-9.32 10-13 



100 



Mo 



present 
ref. [g 
ref. ra 



4.47 10-13 
2.05 10-13 
6.77 10-1^ 



1.93 10 
1.05 10 
3.28 10 



-12 
-12 

-14 



6.58 10-* 
3.50 10-* 
2.91 10-9 



-4.22 10-13 
-1.61 10-13 
-5.64 10-15 



1.32 10- 

6.48 10- 

-1.11 10 



10 

11 
-11 



-2.03 10-12 
7.03 10-13 
2.45 10-1^ 



128 



Te 



present 



ref. [|4 
ref. m 



3.60 10-1^ 
3.36 10-1^ 
1.36 10-1^ 



1.12 10 
7.39 10 
4.32 10 



-14 
-15 
-15 



3.14 10-9 
1.50 10-9 
8.17 10-10 



-1.10 10-1^ 
-4.86 10-15 
-5.24 10-15 



1.42 10-11 

9.46 10-12 

-4.73 10-12 



-1.25 10-1^ 
-1.87 10-15 
-8.51 10-15 



130 



Te 



present 
ref. [P 



ref. [ig 



7.83 10-13 
5.34 10-13 
3.02 10-13 



1.97 10-12 
1.05 10-12 
7.44 10-13 



5.66 10-* 
2.25 10-* 
1.61 10-* 



-5.19 10-13 1.75 10-10 
-2.17 10-13 9.10 10-11 
-2.28 10-13 -6.24 10-11 



-2.34 10-12 
-4.13 10-13 
-1.49 10-12 



The nuclei ^^Ge, ^"^Se, lo^Mo, i2*Te and i30Te have been evaluated within an eleven di- 
mensional model space including all single particle orbitals of oscillator shells 3huj and 4huj 
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plus the Ohg/2 and 0/iii/2 orbitals from the 5huj oscillator shell. In the case of ^^Ca we work 
in a seven dimensional model space including all the orbitals in the major shells 2hu and 3huj. 
Here, the experimental single-particle energies have been used for the orbitals lpi/2, O/5/2, 
lP3/2, 0/7/2, lsi/2 and 0d3/2, while for the remaining orbitals a single-particle energy spacing 
of huj = 41 A~^^^ MeV has been assumed. Finally, both the T = 1 and T = proton-neutron 
interaction strengths in the particle-particle channel have been set by following the recipe in- 



troduced in ref. 28 



Our results for nuclear matrix elements are compared in Table |^ with those obtained by 
MBK and PSVF. In both works configuration spaces similar to ours were employed, and the 
FNS effect was included in the way we have done it (see eq. (p3l)). Yet there are two differences 
that could in principle be important: i) instead of the 6-ioTce, they have used the G-matrix 
(derived from the nucleon-nucleon potential) as the residual interaction, and ii) their correlation 
function is not that given by eq. (^). In spite of these dissimilarities, our results concord 
surprisingly well with those obtained by MBK, except for M^ and Mp. [| The major difference 
is found in ^^^Mo, but we know that this is a "difficult" nucleus from the nuclear structure 
point of view, because of the collapse of the QRPA in the physical region of the particle-particle 
T = strength. Moreover this effect is amplified by the SRC. The agreement with the PSVF 
calculation is only good in the case of the Gamow- Teller moments. Note that the last can also 
be say for the concordance between the MBK and PSVF results. 

The coefficients Cij, defined in eq. (^) and evaluated with the matrix elements given in 
Table ^, are compared in Table |^. Kinematical factors from ref. have been used in our 
calculations. Obviously, all the above mentioned differences between the matrix elements are 
reflected on the calculated Cij values. However, the spread between the entries in the same row 
in Table |^ is smaller than the spread of the values in Table |^. This is because the effect of the 
matrix elements Mj, and Mp is comparatively small. 

Finally, Table |^ gives the constraints on the Majorana neutrino mass and the right-handed 
coupling constants, deduced from the most recent experimental bounds for the /3/5oi/ half-lives, 
and the present evaluation of the nuclear matrix elements. It should be kept in mind that in 
doing so we have used the bare value g^ = 1.254 for the axial-vector coupling constant, and 



^Except the tensor moment Mj,, the MBK matrix elements agree remarkably well with those obtained by 



Tomoda and Faessle r |13 . These authors have not evaluated M^, since a negligible small value for it was 



obtained previously p7| in a projected mean-field approach. 



13 



that the upper hinits for the lepton violating terms shown in Table do not simply scale as g"^. 



Table 7: Experimental half-lives for the neutrinoless double beta decay and upper limits on the 
Majorana neutrino mass (m,^), and the right-handed current coupling strengths (A) and (rj). 



Nucleus To,y{exp) [yr] |(^i/)| [eV] 



l(A)| 



1(^)1 



^'Ca 


> 1.1 10^^ ") 


< 


15 


^^Ge 


> 1.2 10^^ ^) 


< 


0.51 


''Se 


> 2.7 10^2'^'^) 


< 


5.3 


^'>^Mo 


> 5.2 10^2 ^M 


< 


3.4 


i28Te 


> 7.7 10^4 "^^ 


< 


0.97 


ISO^g 


> 8.2 10^1 ^) 


< 


6.4 



-9 



< 1.1 10-5 < 7.9 10-* 

< 8.1 10-^ < 3.2 10 

< 5.7 10-6 < 3.2 10-s 

< 3.2 10-6 < 1.7 10-^ 

< 3.4 10-6 < 6.4 10"^ 

< 7.9 10-6 < 4.6 10"^ 



"•* (laboratory data) ref. pO 
''^ (laboratory data) ref. |^| 



'^^ (laboratory data) ref. [22 
'^^ (geochemical data) ref. [|^ 



(laboratory data) ref. ||2^ 
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5 Concluding remarks 

Nuclear moments for the neutrinoless double beta decay have been evaluated numerically for 
several nuclei, using the formalism that we have recently developed. Simple analytic expressions 
for the j3[3 decay of '^^Ca, that follow from this formalism in the single mode model, are also 
presented. The results shown in Table ^ are useful, not only for testing the full numerical 
calculations, but also for checking the consistency with other formalisms [^ ^ §, P]. In fact, it 
would be highly desirable to find out whether these formalisms lead to numbers shown in Table 
^ This would be a simple and definite test for all nuclear matrix element except for the Mp. 
However, even in this case, a simple model can be framed for confronting different formalisms 
with each other. 

The present work differs from similar QRPA studies, not only in the (3(3f)jj formalism, but 
also in the residual interaction. Namely, we have used a simple 5-force, instead of the G-matrix 
that is currently employed. The fact that our results are equivalent to those obtained by 
MBK and PSVF clearly shows that the (3(3qjj half-lives are not very sensitive to details of the 
nuclear force. In other words, the so called "realistic interactions" are in no way the panacea 
for the nuclear structure evaluation of the j3(3 processes, as it has been proclaimed by some 
authors for a long time. The reason for that is the crucial role played by the restoration of the 
isospin and SU(4) symmetries, produced by the residual interaction, in tailoring the Fermi and 
Gamow- Teller transitions strengths, respectively f^. This restoration mechanism is not limited 
to the RPA-like models p|, ^, ^ ^, but also occurs in the shell model calculation. (One 



simple example has been discussed in Sec. 3.) We are convinced that the nuclear structure 
issue, involved in /3/5-decays will still keep us busy for a long time, and that definitively it 
cannot be solved by a simple minded employment of "good" residual interactions. Neither the 



renormalized nor self consistent RPA methods are able to fix this problem |]32 . 

An alternative technique for examining the (3(3 transitions could be the relativistic RPA, 
which has recently been successfully applied to the description of the isobaric analogue and 



Gamow- Teller resonances in closed shell nuclei [^. We are planning to analyze the conse- 
quences of such an approach. 
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